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jy^ , Abstract. The totally asymmetric simple exclusion process (TASEP) with a 

localized defect is revisited in this article with attention paid to the power spectra 
of the particle occupancy N(t). Intrigued by the oscillatory behaviors in the power 
spectra of an ordinary TASEP in high/low density phase(HD/LD) observed by Adams 



' et al. (2007 Phys. Rev. Lett. 99 020601), we introduce a single slow site with hopping 



CZ3 



rate q < 1 to the system. As the power spectrum contains time-correlation information 
of the particle occupancy of the system, we are particularly interested in how the defect 
affects fluctuation in particle number of the left and right subsystems as well as that of 
the entire system. Exploiting Monte Carlo simulations, we observe the disappearance 
of oscillations when the defect is located at the center of the system. When the defect 
is off center, oscillations are restored. To explore the origin of such phenomenon, we 
use a linearized Langevin equation to calculate the power spectrum for the sublattices 
and the whole lattice. We provide insights into the interactions between the sublattices 
coupled through the defect site for both simulation and analytical results. 
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1. Introduction 

There has been increasing interest in the investigation of complex phenomena, 
most prominently in problems inspired by biology in which a network of reactions 
are coordinated with surprising efficiency and accuracy. At the heart of these 
problems, however, are smaller fundamental ingredients and various pathways. Physics 
is notably well-suited for approaching such complex problems from a "reductionist" 
perspective: Using simplified models to capture the essential constituents and describe 
the fundamental behaviors. One of the paradigmatic systems in non-equilibrium 
statistical mechanics is the totally asymmetric simple exclusion process(TASEP). Using 
a matrix-product formalism[U E] or recursion relations)!], the steady state of TASEP 
with open boundaries can be obtained exactly. TASEP has been applied to model 
both transcription ([U [5]) and translation([6j [7J El |9]) processes in protein synthesis 
as well as molecular motors (pill [IT]), vehicular traffic ([121 02]) and inhomogeneous 
growth processes ([HI [15]). Computationally affordable, TASEP also shines light on 
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properties of complex systems that are experimentally unattainable at the moment. The 
steady state properties of TASEP have been extensively studied both in the abstract 
and with applications. Despite its simple construction, TASEP displays non-trivial 
behaviors from a theoretic point of view such as phase transitions in one-dimension 
with short-range interaction. Recent investigations on TASEP in frequency space 
[161 IZl HH1 UHl E01 [2H [221 [231 121] provide insights into the rich dynamic behaviors of 
this model. Of more relevance to this article, the authors of [22| 121] look into the power 
spectrum of the time-dependent particle occupancy of TASEP and uncovered even more 
surprises when oscillations are observed for TASEP in high/low density phases. 

In this article, we focus on the power spectrum of TASEP with local inhomogeneity, 
especially the interaction between the two sublattices separated by the inhomogeneity. 
We outline the biological motivation in Section [LT] and review previous works pertaining 
to this study in Section 11.21 We define our model and provide simulation results in 
Section [2j Using Langevin dynamics, we present our theoretical results which provide 
insights into the interactions between the sublattices in Section [3] We conclude in 
Section @] 

1.1. Biological Motivation 

In a living organism|25], proteins are not only the building blocks of cells (comprising 
most of a cell's dry mass), but they perform nearly all cellular functions as well. To 
optimally synthesize proteins, namely first transcribing the genetic information from 
DNA into a messenger RNA (mRNA) and later translating mRNA into a functional 
protein, is of crucial value in vivo for the organism and in vitro for pharmaceuticals. 
We devote our attention to the translation process here. 

The template of translation, mRNA, is composed of nucleotide triplets called 
"codons" that contain information for the desired amino acids. As there are four distinct 
nucleotides(A, U, C, G), there are 64 codons. Except for three stop codons(UAA, UGA, 
and UAG) that signal the termination of translation, the rest "code" for one of the 
twenty amino acids with degeneracy ranging from 1 (e.g. AUG for methionine and 
UGG for tryptophan) to 6 (e.g. CGN and AGR for Arginine). In prokaryotes such 
as a bacterium, translation initiates when ribosomes - large molecules acting as an 
assembling machinery - encounter the start codon AUG. During the elongation stage, 
ribosomes move along the mRNA and incorporate transfer RNAs(tRNA) with correct 
amino acids attached. Signaled by one of the three stop codons, translation terminates 
and the amino acids are later joined by peptide bonds to form different sequences and, 
with proper folding, result in various final protein products. 

Codons coding for the same amino acid are termed "synonymous." Even among 
synonymous codons, different tRNAs may be employed. Therefore all codons are 
translated, in principle, at a different rate depending on the availability of the cognate 
tRNA's. When a certain tRNA becomes scarce or even unavailable, the ribosome awaits 
on the corresponding codon ( "slow codon" ). This slows down the translation process and 
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can lead to translation error or frame-shifting which potentially results in an erroneous 
protein product or terminates translation prematurely. In light of the tRNA cellular 
concentrations in Escherichia coli [26J which can differ by nearly an order of magnitude 
even among synonymous codons, we look into the effect of having one slow codon 
at different locations along the mRNA. Exploring the power spectrum of such modied 
TASEPs is expected to provide information on the fluctuation in the timescale of particle 
transport through a finite system, a feature of biological relevance. 

1.2. Review of previous works 

TASEP comprises of a one-dimensional (ID) lattice of L sites and particles of size 
1 in units of lattice sites moving uni-directionally with hard-core exclusion. The steady 
state properties of an ordinary TASEP with open boundaries, with entry/exit rates a/ [3 
and particle hopping rate 7 = 1, have been extensively studied([UEll3])- Given different 
a and f3, the system can be found in a low density (a < [3, a < 1/2, "LD"), high density 
(p < a,p < 1/2, "HD"), maximal current (a > 1/2,(3 > 1/2, "MC") and a shock phase 
{a = (3 < 1/2,"SP"). 

It is not until recently that the dynamics of TASEP are getting more attention [161 
HZl HH1 UHl EQl EH [221 EH [24]. In [Hj], a combination of domain wall theory and 
Boltzmann-Langevin theory well characterized the dynamics of the LD and HD phases. 
Soon after, the complete power spectrum of the transition matrix of a partially 
asymmetric exclusion process was obtained using the Bethe ansatz[20l EI], providing 
an accurate description of the short-range particle-hole correlations. The authors of [22] 
investigated the effects of system sizes contained in the power spectra of the time series 
of system occupancy N(t), which provides insight into the time-correlation of N(t) and 
how long it takes the noise to propagate through out the entire system. Defining a 
Fourier transform N(u) = Ylt=i N(t)e tuJt where u = 2irm/T, T is the length of each run 
and m G [0, T], the authors obtained the power spectrum by averaging over 100 such 
realizations: 

/M = (|iVM| 2 >. (i) 

When looking at the HD and LD phases, they found marked oscillations of 
that damp into power laws. They further explored the origin of such behavior by taking 
the continuum limit and using a stochastic equation of motion. Their analyses and 
simulation results reveal the oscillatory minima capture the system size and are located 
at n ■ 2nv/'h (or n ■ vT/h in m) where n is an integer, v = 1 — 2p and p is the average 
bulk density. Note, the quantity u/L is simply the time the fluctuation takes to traverse 
the lattice. 

Their results are useful in that many physical systems, e.g. mRNA with several 
hundred codons, are far from the thermodynamic limit. The timescale at which noises 
traverse the entire system can be significant for particle transport. In the case of cellular 
protein synthesis where multiple mRNA's are translated simultaneously, knowing the 
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finite size effect on the fiuctuations of total occupancy could shine some light on how 
mRNA competes for the "translation machinery," i.e. the ribosomes. 

Our study extends the investigations in [22] to studying the power spectra of TASEP 
with one localized defect (or "slow codon" in biological sense) where the hopping rate 
q < 7 = 1. Since an mRNA consists of a series of completely inhomogeneous 7's, 
studying the effect of one defect to the overall particle throughput serves as an ideal 
starting point. We explore both the entire system and the two subsystems joined by 
the defect. Our simulations show that when the defect is off the center of the lattice, 
the power spectrum of the entire system reflects such "bias" while each sublattice can 
still be characterized by modifying the theory for an ordinary TASEP. When located 
in the middle of the lattice, the cross correlations between the sublattices eliminate 
the oscillations, leaving a power law decay. The continuum theory in [22] cannot be 
readily modified to account for this case, so we opt to use discrete space and time units 
to analytically explore the power spectra for all locations of the slow site. Our results 
furthered the work by Adams et al. [22] in the following sense: a) Since the defect can be 
located anywhere throughout the lattice, the left and right subsystems can have different 
sizes and therefore a different timescale for a hole-particle pair to propagate through. We 
provide both simulation and analytical results to investigate such difference; b) The left 
and right subsystems are coupled through the defect, leaving this site and its rightmost 
neighbor crucial in theoretical treatment. More careful considerations are needed which 
we present in the later section of this article. 

2. Model and simulation results 

We define our system on a ID lattice of L sites with open boundaries. L = 1000 

for most simulation data unless otherwise specified. Particles, injected from the left 

end with rate a and taken out at the right end with /3, move along the lattice with 

7 = 1 except at the k th site where a defect is introduced with hopping rate q < 7. 

Using a mean-field approximation |27J, one finds when the defect site is far from the 

system boundaries, the left and right sublattices can be viewed as two TASEP's joined 

respectively by exit and entry rate of 

q 

q eS = - ~ q when g<Cl. (2) 

To ensure a flat density profile without loss of generality, a and (3 are chosen to equal 
(1 ~q e «)- 

We use a random sequential updating scheme and define one Monte Carlo step 
(MCS) as making L + 1 attempts to update the lattice, giving on average one chance 
to update each lattice site as well as introduce a new particle. All measurements are 
taken after 10K MCS to ensure the system has reached steady state. N(t) is collected 
every 10 MCS for t = 1,2,...T where T = 10 4 (10 5 MCS) for most of the data. We 
adopt the definition of power spectrum in Eq.flTJ, but with a 1/10 5 normalization. The 
results are averaged over 100 realizations. When investigating the effect on the power 
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spectrum due to one defect at different positions along the lattice, we measure the power 
spectra for both the entire system Itot(w) and the left and right subsystems separated 
by the defect, Iz(u>) and Ir(lo). As in previous studies, oscillations are found in the 
LD/HD phases. However unexpected and intriguing results emerge when the location 
(k) and the "strength" (q) of the defect are varied. We shall elaborate these results in 
the following sections. 

2.1. Defect off center: k ^ L/2 




200 400 600 J00 1000 1200 1400 



Figure 1. Power spectra of the entire lattice (top, black online), left (middle, blue) and 
right (bottom, red) sublattices with L = 1000, k = 333 and q = 0.001. w = litmjT. 



Let us first introduce a single defect of hopping rate q at site k ^ L/2. The 
average bulk density of the system becomes: 

k L k . . 

PTOT = jPL + ~^—pR- (3) 

When the defect is off-center, Ptot is found in either HD or LD phase, while the 
left sublattice is in HD with p L = \ — q cS and right sublattice in LD with p R = q cS . Once 
a hole-particle pair is created at site k, it introduces fluctuations in the hole density to 
the left and particle density to the right. As the defect is off center, the fluctuations in 
the left propagate with vl = 1 — 2g efi = v for a distance of k and the fluctuations in the 
right with vr = v for a distance of (L — k). We therefore expect the minima in Il{oj) 
and Ir(cj) to be at n-vT/k and n-vT / (L — k), respectively, since each sublattice may be 
treated as an independent TASEP as in [22]. Concerning the entire lattice, on the other 
hand, the power spectrum measures the difference in timescale between the hole and 
particle fluctuations propagating through the left and right subsystems. The particle 



Power spectra of TASEPs with a localized slow site 



6 



and hole fluctuations propagating through the sublattices do not cause a fluctuation 
in the overall particle density; only when one of the fluctuations remains following its 
counterpart leaving the system does the total lattice density change. Therefore, the 
minima of Itot{^) are located at n ■ vT / (L^ — Lr) = n ■ vT / (L — 2k). 

FigJU displays the simulation result when k = L/3 = 333 and q = 0.001. It is 
reasonable to take q eSl ~ q — 0.001, yielding v = 0.998. As expected, both Il{oj) and 
Ir(oj) demonstrate oscillations of which the minima well characterize the subsystem 
sizes: 

Minima of left sublattice: = n ■ (0.998) 10 5 /333 w n ■ 300 
Minima of right sublattice: = n ■ (0.998)10 5 /666 w n ■ 150 

We also observe oscillations in Itot{u>). In this case, \L L — h R \ = h L and hence we 
expect the minima in Itot{u) to be at the same location as those of Il{w), which is 
precisely captured in Fig. [TJ 

2.2. Intermediate range of uj 




500 1000 5000 
m 

Figure 2. I T OT(m) with m G [500,4000] for L = 1000, k = 333 and q = 0.001. 
to = 2irm/T. 

In addition to the previous simulation findings, we want to comment on, without 
going into much detail, the behavior of Itot{<^) demonstrated in the intermediate range 
of u, i.e. 1 « m < T when oscillations are present. Consistent with the simulation 
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and theory results in [22] , here we retrieved the power law decay in the power spectrum 
of the entire system. For a system size of 1000, an exponent of -1.5 decently captures 
the decay (Figj2]). Upon closer examination, the rate of decay varies with the system 
size. As discussed in [22J, the exponent approaches -2 if the system is large enough 
to allow the oscillations to damp completely. Our simulation confirms this conclusion. 
The continuum approach adopted in [22], on the other hand, is not suitable to analyze 
the system with a defect. Our theoretical treatment turns to the discrete space and 
time units for the low-frequency regime. We leave the intermediate u regime for future 
pursuit. 

2.3. Defect at center :k = L/2 



NO oscillation 




200 400 600 ^00 1000 1200 1400 



Figure 3. Power spectra of the entire lattice (top, black online), left (middle, blue) and 
right (bottom, red) sublattices with L = 1000, k = 500 and q = 0.001. u = 2nm/T. 

Now we turn to a special case when the defect is located at the center of 
the system (A; = L/2). Due to particle- hole symmetry, the system is always half-filled 
(Ptot = 1/2)- The left and right sublattices again are in HD and LD, respectively and 
the minima in the power spectra are at the same location:?-;, • vT/k. When looking at 
Itot{^)i what behavior shall one expect when — Lr = 0? Now that the hole-particle 
pair is created right in the middle, it takes on average the same time for the hole and 
particle to exit the system. However, it is the fluctuation of the timescale that Itot{^) 
tries to capture. 
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Through our simulation data in FigfJ] where k = L/2 = 500 and q = 0.001, we 
again retrieve the oscillations in and Ir(oo) of which the minima are both at: 

n ■ (0.998)10 5 /500 « n ■ 200 

However no oscillation is observed in J TO t(^)- The origin of such behavior is explored 
in Section [3j 



3. Theory 

In this section, we present an explanation for the previously discussed simulation 
results. Our approach utilizes a Langevin equation for the site density p(x, t), which we 
linearize about its average p(x) [22], [24]. In addition, the theory uses discrete space and 
time to naturally include UV cutoffs, which the authors of [21] previously used for the 
constrained TASEP. Since mRNA's typically contain several hundred codons, we avoid 
taking the continuum limit that was used previously in [22j. Further, we do not take 
the Fourier transform of the spatial component as in [22] and [23] due to the breaking 
of translational invariance by the slow site. Instead, we solve a set of L linear equations 
obtained from the Langevin equation for each site. We also choose the timescale to be 
a single update attempt with T being the number of updates in a run as in [23]. This 
choice is in contrast with what is chosen for the simulation, where data is taken every 
10 MCS. To account for the difference, we must sum over the harmonics [23] of our 
theoretical results to make a meaningful comparison. 

Using a similar approach as [24J, we begin with the Langevin equation for the site 
density: 

d t p(x, t) = P(^-l^)[l-P(^^-P(^*) [1-^ + 1.*)] + ax _ M) _ t) (4) 

where x ^ {1, k, k + 1,L} and time is on the order of an update. For those sites, we 
have: 

MM) ^[i-P(M)]-P(M)[i-P(2. t )] +;(0t) _ ;(M) (5) 
g , k () = Pit ~ h ''I 1 Z P(M)1 - qpjk, t)[l - p(k + M)] 

+ S ( k-i,t)-£{k,t) 

d tP (k 1 1;t) _ gp( fc ^)[ 1 -p( fc + 1 ^)]-p( fc + 1 ^)[ 1 -p( fc + 2 ^)] (7) 

+ Z(k,t)-{;(k + i,t) 

dtp{h t) = p(L - 1 , t) [1 - p(L, gj - (3p(h, t) [1 - p(L, t)] (g) 

+e(L-i,t)-e(L,t) 

where q is the hopping rate for the slow site and £ is assumed to be homogeneous, 
uncorrelated white noise associated with hopping from site x to x + 1 with zero mean 
and variance A such that t)^(x', t')) = A5 XyX i5 t ^ [24J. This assumption for the noise 
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will later be validated by the agreement of our analytical result with the simulation 
results. The 1/(L + 1) on the right-hand side is due to the probability of choosing a 
particular site to update. 

To simplify the calculation, we choose a and /3 such that the profile is flat (i. e. no 
x dependence) for the two regions separated by the slow site. This is also consistent 
with the choice of parameters in the simulations. Thus, we have 

p(x) = p L xe [1, k}; pi = = a; (9) 

p(x)=p R xe[k + l,L}; p R = -^— = f3 (10) 

Note, pl + Pr = 1 with this choice of parameters. Particle/hole symmetry emerges with 
this choice of parameters as Pl + Pr = 1- We then define the fluctuations of p(x,t) 
about its average to be 

(p(x,t) = p(x,t) - p L xe[l,k] (11) 
<p(x,t) =p(x,t) -p R xe[k + l,L] (12) 

For small values of q, the particle spacing on the right sublattice is so large that 
the interactions do not contribute significantly. Similarly on the left sublattice by 
particle/hole symmetry, the hole spacing should be large enough that they do not 
interact with each other. These interactions appear as quadratic terms in tp(x, t), which 
we will assume to be small and dropped in the following. 

At this point, we will treat the left and right sublattices separately. Rewriting the 
Langevin equations and keeping only terms linear in tp, we have 

MMi ^y' tflM-flM) (i3) 

MM) - - H)MX - M) Z^M+Mit+IA +{(,-!,«) -{(,,,) (14) 
x e [2, k - l] 

m*. o = (1 = h)ip(k = M) I - + PMk + ~ +(( k ~ m) -((k,t) (is) 

M* + 1,0= (i^MM^glM+jj^+M + {(tj () _ {(Jb + lf () (16) 

wm) - (1 = feV(x = u) L f " + ^ + M> + - m> - «*. o 

ie[fe + 2,L-i] 

WM ) = ^M^ + ; (l _ M) _ {(Li() (1 8) 

Taking the Fourier transform to be 

V (x,t) = J2e- lu ' t ^(x,co) (19) 

^,^) = ^e^V0M) (20) 
t 
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u= 2 -^ me[l,T] (21) 



we arrive at 



C(u#(l, u) - |^;0(2, w) = r/(l, w) (22) 
^ L -(p{x - 1, w) + C(w)y5(x, w) - — ^—ip(x + l,o;) = ^(x, u) (23) 



L+ 1 rv ' y * w-v > / L + 1 

x 6 [2, Jfe - 1] 

' ~ ^ L -(p(k -l,u) + C(cj)<p(k, u) - F^r^fc + 1, w) = fj(k, u) (24) 



L+ 1 rv ' ' >v /n ' ' L+ 1 

1 ~ w) + CM£(ife + 1, w) - T^—V{k + 2,u) = fj(k + 1, w) (25) 



L+l v ' L+l 

x 6 [Jfe + 2,L-1] 

1 — fl» 

-0(L- l,w) + CM<^(L,w) = fj{L,u) (27) 



L + l 



where 



fj(x,u) = i(x - l,w) - (29) 
or written more compactly, 

M<^ = 77 (30) 
where y5 and fj are vectors in x. The LxL tridiagonal matrix M has entries: 

M u = C(w) (31) 

Mii+1 = "LTT ^t 1 '* -1 ] ( 32 ) 

M ii+1 = - s ^ I ie[k,L-l] (33) 

M ii-i = -]f-vY ^ G t 2 ' A; + 1 ] ( 34 ) 
Mll ~ 1 = ~i7+i ie \ k + 2 ^ v \ ( 35 ) 

where we have made use of Pl + Pr = 1- By inverting the matrix (M _1 = §), we have 
the <£>'s in terms of the 77's, namely 

L 

<p(x, u) = ®xvV(y, w) (36) 
The inverse of this type of tridiagonal matrix is given in [28J. 
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We are interested in power spectrum of I = (|iV| 2 ), so first we need N. N is related 
to (p by summing over the spatial coordinate. Thus from (122]) . we have 



N TOT (u) 



PL 

L+ 



i [£(1, w) + w)] + i(h, co) - |(0, w) 



1 - e 1 ' 



lXi + 1, w) - p fl £(fc, w) - w)] + |(0, w) - 



lTt w) - Pi?y?(& + 1, w) - Pl<p(L, w)] + w) - £(L, w) 



(37) 
(38) 
(39) 



for the entire lattice and the left and right sublattices, respectively. Upon squaring 
N and averaging over the noise to calculate the power spectra, we still need to sum 
over the harmonics to make meaningful comparisons with the simulation data [21]. The 
harmonic sum is necessary since the simulation data is not recorded at every update. If 
we take data every £ updates, then 

i-\ 

I{m) = Y,I^m, z ) (40) 



where 



z=0 

2tt 



m + 



zT 



(41) 



Finally, the only remaining parameter in the theory is the variance of the noise A. 
Instead of leaving it as a fit parameter, A is obtained from the microscopic description 

m 



A (L + l)(l + g) 2 
which is simply the current divided by L + 1. Thus, the final result is 

e-i 



(42) 



'TOT 



(m) = ^ 



q 



- (L + l)(l + g) 2 T(2-2cosu; miZ ) 



(43) 



PL 



L 



-j {(tp(l,oj mtZ )0*(l,u mtZ )} + ((p(L,u mtZ )(f*(L,u m:Z )) 



L + l 

I r(m) = > — — 

M ; ^ L+l 1 



+ 23^(£(l,w m ,,)£*(L,w m , z ))]} 



i-i 



$l{([0(l,uj m , z ) + ^,eJm, z )] r(L,w m , 2 ) -r(0,w m , 2 ) }\+2 

(44) 



2=0 



+ l)(l + g) 2 T(2-2cosw m , 2 ) 



x 



PR 



L + 



- J {(ip(k, U m:Z )(p*(k, U m ,z)} + (<f(k + l,UJ m)Z )<p*(k + l,UJ m , z )) 



- 25K [(<p(k, u m)Z )ip*(k + 1, U} m ,z))]} 



PL 

L+l 



{(p{l,u m ^)(p*{l,u miZ )) 
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+ 
+ 
+ 

J «M = E 7l 



ZpLpR 



L + l) 2 
L + 
L + l 



5i {([</?(&, w m , z ) - (p(k + l,w m ,2f)] 0*(l,^m,«))} 



>} 



+ 2 



2=0 



+ l)(l + g) 2 T(2-2cosw mi2 



(45) 



x 



PR 



L 



- 2K [(<p(fc, uj m , z )(p*(k + 1, w m , z )}]} 
2pLPi? 



Pi 



L+l 



(<p(L, w mi2 )v?*(L, o; m , 2 )) 



L + l) 2 

2p_i? 
+ 

2p L 
L + l 



3? {([<£(& + 1, w m , 2 ) - w^)] (f*(L, w m ,g))} 

ft\(v(L,cu m , z ) \£*(L,LJ m<z ) -£*(k,U} m , z j\)\ + 2 



l>} 



where 



(47) 
(48) 
(49) 



(0(a,u)<p*(b,u)) =2j2®a X % x ~ [ § -%-D + S a (*-i)Sy (46) 

x=l x=2 

(<p(a,u)£*(0,u)) = S„i 
(<p(a, w)|'*(L, w)) = -S aL 
(<p(a, w)|*(6, w)) = S a (6+i) - Sa6; 6 7^ 0, L 

and the dependence of S on u has been suppressed. 

Before making a comparison with the simulation results, we will focus on the special 
case of k = L/2. In this case, both sublattices have the same lengths, while we still have 
the particle/hole symmetry previously mentioned. Using these symmetries, we reduce 
the number of unique elements in the § matrix. Specifically, we have 

= §(L-;+i)(L-j+i) (50) 

Further, the cross correlation between the first and last site {<p(l,u)<pD, which 
only appears in the power spectrum of the entire lattice, is a real quantity. The cross 
correlations in the sublattices' power spectra are not real. In the previously studied 
homogeneous TASEP, the oscillations appeared mathematically as the real part of a 
complex function [22]. Here, it is less clear how the oscillations emerge, as opposed to 
the simpler expression in [22J. However, we speculate that the real correlation terms 
are the source of the disappearance of the oscillation, and we are currently investigating 
this idea. 
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Using equations (pE3"j) . (03J, and (HS]), we compare the theoretical result to the 
simulation data. For low q values, the agreement is remarkably good. In figure HI we 
show the comparison for L = 1000, k = 333, and q = 0.001. As shown in the figure, the 
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Figure 4. Comparison of simulation and theoretical power spectra for L = 1000, 
k = 333, and q = 0.001. 

oscillations for all three power spectra are correctly predicted for the entire range of m. 
The theory also predicts the vanishing of the oscillations for the total power spectrum, 
as shown in figure EJ for L = 1000, k = 500, and q = 0.001. Again, the theory correctly 
predicts the oscillations of the two sublattices with very good agreement. Therefore, 
the theory incorporates all the necessary information needed to accurately reproduce 
the power spectra for small q values. 

As we increase the q value, a difference between the simulation data and the 
theoretical result begins to appear. This difference is shown in figure [6] for L = 1000, 
k = 333, and q = 0.2. While the theory does correctly predict the m value of the 
minima, the damping of the oscillations is greater in the simulation than in the theory. 
Physically, the particles/holes begin to interact more frequently (via the exclusion) at 
larger q values. However, we employ a linear theory void of these interactions. The 
effect of the interactions is an enhancement of the effect diffusion coefficient [22j 129] . 
which controls the damping of the oscillations [22] 129] . Therefore, it is no surprise that 
the theory is not in perfect agreement with the simulation data. Since the interactions 
terms (terms quadratic in ip) are assumed to be small, perturbation theory should be 
utilized to account for the difference, not only here but for the homogeneous TASEP 
[22] and other modified TASEP models [21] as well. 
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Figure 5. Comparison of simulation and theoretical power spectra for L = 1000, 
k = 500, and q = 0.001. 
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Figure 6. Comparison of simulation and theoretical power spectra for L = 1000, 
k = 333, and q = 0.2. 



4. Summary and outlook 

In this paper, we explored the power spectrum associated with the total occupancy 
of a TASEP with a single slow site. The slow site (at position k with rate q) served as 
a bottleneck for particles passing through from the left to the right sublattice. Using 
Monte Carlo simulations, we observed oscillations, similar to those seen in the HD /LD 
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phases of the homogeneous TASEP [22], for the two sublattices separated by the slow 
site and for the entire lattice. However, the oscillations for the total power spectrum 
vanished when the slow site was place in the center of the lattice, while the oscillations 
remained for the sublattices. A theoretical approach utilizing a linear Langevin equation 
(with discrete space and time) for the density correctly captured this disappearance 
for k = L/2 as well as the oscillations for the other cases. The theory had excellent 
agreement when compared to the simulation results when the hopping rate for the slow 
site q was small. As the rate was increased, the overall agreement was not as impressive, 
but the locations of the minima of the oscillations for the simulation and theory were 
the same. 

While the theoretical result presented is not very transparent, we can glean some 
insight into the disappearance of the oscillations for the power spectrum of the entire 
system when k = L/2. Specifically, we found the cross correlation between the density 
fluctuations of the first and last site is a real quantity, as opposed to having an imaginary 
component as well. For the homogeneous TASEP [22], the oscillations emerged from 
the real part of a complex function which we lack in this special case for the single slow 
site in the center. We speculate that the cross correlation being a real function signals 
the disappearance of the oscillations. 

Open questions remain beyond this study. One immediate question is the 
connection between our theoretical approach and the one in [22]. While there is good 
agreement between our theory and simulation results, the theory lacks the simple, 
intuitive explanation demonstrated in [22]. By taking the appropriate limits, we would 
like to have a similar simple expression for the oscillations that can be used to show the L 
and k dependence of the minima for the power spectrum of the entire system, especially 
in the k = L/2 case. Finding such a connection would not only expand our understand 
for a slow site, but other extensions [301 EI] as well. Another question pertains to the 
power law decay found in the k = L/2 case power spectrum for the entire system. The 
simulation data shows a u~ 3 ^ 2 decay after some crossover scale. We wish to explore 
the L dependence of the crossover as well as explore the physics behind the power law. 
Lastly, what is the effect of two or more slow sites [23 E] arid extended objects [31], [32] 
on the power spectrum? 

Many biological extensions of this model exist that would more closely realize what 
happens in nature. For the translation process, these extensions include coupling a 
single slow site with a finite pool of particles [33] to more accurately model the finite 
resources in the cell and multiple TASEPs [30] to capture the effects of competition 
for resources in a cell. Also, one should have completely inhomogeneous hopping rates 
to model the various codon elongation rates. By studying the individual elements and 
their coupled effects, we gain a better understanding of the translation process in protein 
synthesis. Outside of the context of protein synthesis, many other biological systems 
may be studied through the perspective of physics. From social networks to predator- 
prey dynamics, much work is left to be done to understand nature around us. 
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